Quantum Renormalization of the Spin Hall Effect 
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By quantum Monte Carlo simulation of a realistic multiorbital Anderson impurity model, we study 
the spin-orbit interaction (SOI) of an Fe impurity in An host metal. We show, for the first time, 
that the SOI is strongly renormalized by the quantum spin fluctuation. Based on this mechanism, 
we can explain why the gigantic spin Hall effect in Au with Fe impurities was observed in recent 
experiment, while it is not visible in the anomalous Hall effect. In addition, we show that the SOI 
is strongly renormalized by the Coulomb correlation U. Based on this picture, we can explain past 
discrepancies in the calculated orbital angular momenta for an Fe impurity in an Au host. 

PACS numbers: 71.70.Ej, 75.30.Kz, 75.40.Mg 
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When a magnetic impurity with d orbitals is put into 
a normal metal with s- or p-conduction bands, it forms 
a virtual bound state hybridized with the conduction 
states, leading to a broadening with width A (hybridiza- 
tion energy) which is typically of the order of cV. This 
energy scale competes with the Coulomb energy U be- 
tween the electrons in the d orbitals, which is also of the 
order of eV. The latter tends to produce the spin mo- 
ment, while the former induces quantum fluctuation of 
that spin moment, leading to the spin singlet. The com- 
petition between these two interactions defines the time 
scale, or equivalently the energy scale, of the quantum 
spin fluctuation, i.e., the Kondo temperature Tk, which 
can be much lower than both A and C/ [Ij. In the lan- 
guage of the renormalization group, this is described by 
the energy dependent scahng of the various operators. 
This means that even a weak interaction can be ampli- 
fied in the low energy or temperature scale and compete 
with the much larger energy scale due to the correlation 
U. 

In real systems, the d orbitals have fivefold degeneracy 
and these orbital degrees of freedom have often been ne- 
glected in the analysis of the experimental results. The 
reasoning is that either the crystal field splitting Ae is 
much smaller than the hybridization energy A, or that it 
is in the limit with Ae much larger than the Kondo tem- 
perature Tk. However, the naive comparison between 
the bare interaction strengths is dangerous since these 
are scale-dependent running coupling constants. Similar 
nontrivial behavior can also be expected for the relativis- 
tic spin-orbit interaction. 

Orbital degrees of freedom in the impurity scatter- 
ing lead to intriguing phenomena such as the anoma- 
lous HaU effect (AHE) and the spin Hall effect (SHE). 
A charge current perpendicular to the applied electric 



field is produced in ferromagnetic metals (AHE), while a 
spin current rather than a charge current, is induced in 
senuconductors and metals without magnetism (SHE). 
An extrinsic mechanism of these two effects arises from 
the skew scattering, i.e., spin-dependent deflection of the 
scattered electrons due to the spin-orbit interaction. It 
has been shown in the case of the AHE that the resonant 
skew scattering due to the virtual bound d states leads 
to a large Hall angle, of the order of 0.01 Q], compared 
with the typical value of 10"'^. In the AHE, the spin fluc- 
tuation is quenched, owing to the ferromagnetism or the 
external magnetic fleld, and the ratio of the spin-orbit 
interaction A and A basically determines the Hall angle. 
On the other hand, in the SHE, the spin fluctuation is 
active and the Kondo physics can be relevant to the res- 
onant skew scattering. Therefore, a crucial question is 
whether the Kondo effect and quantum spin fluctuation 
can produce an even larger spin Hall angle compared with 
the AHE. 

An important clue to this question comes from an ex- 
periment by Seki et al. [3| on Au/FePt, in which a spin 
Hall angle of 0.114 was observed. Motivated by this ex- 
periment, some of the present authors studied the Fe im- 
purity in Au by a first-principles calculation [j]. The 
Kondo effect of Fe impurities is a historic problem, with 
a low Tk around 0.4K [5[ , high electric resistance at room 
temperature [l| , and an AHE with Hall angle of the order 
of 0.01 [2|- However, a simple fivefold degeneracy of the 
orbitals has been assumed to analyze the experiments. In 
Ref. Q, on the other hand, the orbital-dependent Kondo 
effect of Fe in Au was proposed to explain the nature of 
the experimentally observed giant spin Hall signals (j]: 
It was argued that the Cg orbitals of Fe are in the Kondo 
limit and t2g orbitals are in the mixed-valence region. 
The enhancement of the spin-orbit interaction by electron 



correlation in the t2g orbitals leads to the giant spin Hall 
effect. However this proposal has been challenged theo- 
retically by Ref . [a] , which suggests an effective 3-channel 
Kondo model, involving local and band electrons of t2g 
symmetry, and also experimentally by the x-ray MCD [7| , 
which obtained a rather small value of the orbital angu- 
lar momentum in contrast to the large value calculated 
inRcf.[4j. 

In order to resolve this confusing situation and quan- 
tify the mechanism for the enhanced spin Hall effect, it 
is essential to treat the quantum fluctuations of the spins 
and orbitals systematically. This is impossible in the 
first-principles calculation, which assumes ordered spin 
and orbital moments. In this Letter, we overcome this 
difficulty by using the Hirsch-Fye quantum Monte Carlo 
(QMC) simulation [8|, combined with the density func- 
tional theory (DFT) [^, U^ , to study the renormalization 
due to correlation effects. First, a single-impurity mul- 
tiorbital Anderson model [ll[ is formulated within the 
DFT for determining the host band structure, the impu- 
rity levels, and the impurity-host hybridization. Second, 
the magnetic behaviors of the Anderson impurity at finite 
temperatures are calculated by QMC. 

The single-impurity multiorbital Anderson model is 
defined as 

-ff = ^ eQ(k)4^^CkQcr + ^ {V^kadl^Cl^aa + H.C.) 
k,a,cr k,Q,^,f7 

?,'^ ? 

u' J 

where cj^Qg. (ckao-) is the creation (annihilation) operator 
of the conduction electron with wavevector k and spin u 
in the band a, dt id^^) is the creation (annihilation) op- 
erator of the localized electron at the impurity site with 
orbital ^ and spin cr, and n^o- = dl^d^a- The host en- 
ergy band eQ(k), the impurity energy levels e^, and the 
impurity- host hybridization Vjka, as the one-body prob- 
lems, can be properly obtained by the DFT calculations. 
U {U') is the on-site Coulomb repulsion within (between) 
the orbitals of the impurity, and J is the Hund coupling 
between the orbitals of the impurity. These many-body 
interactions can be exactly treated by the QMC calcula- 
tions, and thus it becomes possible to accurately study 
the quantum fluctuations of the spins and orbitals of the 
impurity. Considering the parameters used in the pre- 
vious calculations for Fc in Au [j|, and the relationship 
t/ = t/' -h 2J [ill, in our following QMC calculations, 
we use the values of [/ = 5 eV, J = 0.9 eV, and U' = 
3.2 eV for most cases, but we shall vary the values for a 
few cases in order to clarify the role of correlations and 
reconcile with past calculations. 

Our DFT calculations are done by the code Quantum- 
ESPRESSO [H. To calculate the impurity-host hy- 
bridization, we consider the supercell Au26Fc, where the 
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FIG. 1. (color online) (a) Hybridization between the ^ or- 
bitals of an Fe impurity and an Au host. We show the LDA re- 
sult including five d orbitals. (b) Temperature dependence of 
the square of the magnetic moment ((Af|)^), (c) the tempera- 
ture times susceptibility T^j , and (d) the occupation number 
(n^) of the ^ orbitals of the Fe impurity, (b)-(d) are the QMC 
results including three d orbitals. See text for details. 



exchange-correlation interactions are described by the 
Perdew-Zunger local density approximation (LDA), and 
the electron-ion interactions are represented by the Rabe- 
Rappe-Kaxiras-Joannopoulos ultrasoft pseudopotentials. 
4x4x4 Monshorst-Pack k points are used. Kinetic en- 
ergy cutoff for wavefunctions and charge density are 
taken as 30 and 300 Ry, respectively. 

Figure [TJa) shows the hybridization between ^ orbitals 
of Fe impurity and Au host, where the hybridization 
matrix element has the form of V^ka = (v'cl-f^ol^a(k)) 
= ■^Ep,„e''''"aap(k)((p4|7Jo|v'p(n)). Hq is the one- 
particle part of Eq. dD), ip^ is the ^ state of Fe impu- 
rity, and 'I'Q,(k) is the Au host state with wavevector 
k and band index a, which is expanded by atomic or- 
bitals iy9p(n) with orbital index p and site index n, actp(k) 
are expansion coefficients, {ip^\HQ\ipp{n)) are mixing inte- 
grals, and A'' is the number of lattice sites. It is observed 
that, at the F point (k=0), the hybridization value of the 
^ = eg(z^,a;^-y^) orbitals of the Fe impurity is smaller 
than that of the ^ = i23 {^z,yz,xy) orbitals. 

Based on the above DFT/LDA calculation, we can de- 
termine approximate impurity levels e^= -1.9 eV for ^ 
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and 6^= -1.8 eV for S, ~ t^g {xz,yz,xy) 



with zero Fermi energy, where the Coulomb terms in- 
cluded in the LDA calculations have been subtracted as 
in Refs. [li,[il]. The crystal field splitting Ae = 0.1 eV 
is in agreement with the previous LDA calculations [J, Q . 
To study the renormalization due to the correlation 
effect for (i) the crystal field splitting and (ii) the rela- 
tivistic spin-orbit interaction, the magnetic behaviors of 
Eq. ([T|) are calculated by the Hirsch-Fye QMC simula- 
tion for the infinite host system Q. Owing to calcula- 
tion constraints, we simplify it to a three-orbital model 



with one Cg and two t2g orbitals, which captures the 
essential physics. The Cg orbital is arbitrarily chosen 
as z^, and the two t2g orbitals are chosen as pi and 
p_i, where the notation corresponds to the transforma- 
tional properties of <2g orbitals equivalent to effective 
p orbitals ^M,- Pi = ^Tf'-'^'^ ~ W^)^ Pa = ~ixy, and 
P-i = — 75 (^-^ + W^)- The following results are ob- 
tained with more than 10^ Monte Carlo sweeps and the 
Matsubara time step At — 0.25. 

Figures [ll[b)-(d) show the temperature dependence of 
the square of the magnetic moment ((M|)^), the tem- 
perature times susceptibility Tx^, and the occupation 
number (n^), which are defined as Af| = n^^ — n^j^, xe, 

= Jg dr(M|(r)M|(0)), n^ = n^^ + ri^j,, respectively, 
and the impurity levels e^— -1.9 {^~z^) and -1.8 eV 
{^=Pi, P-i) are used. It is found that the eg{t2g) orbital 
has a larger (smaller) magnetic moment, a much larger 
(smaller) susceptibility with a much smaller (larger) 
Kondo temperature, and a larger (smaller) occupation 
number. In addition, if we repeat the QMC calcula- 
tion with the degenerate impurity levels by hand, e^~ 
-1.85 eV for S,=^z^, pi and p-i, nearly the same behaviors 
are observed. Thus it is clear that the orbital-dependent 
Kondo effect comes mainly from the renormalization of 
the impurity-host hybridization by correlations. 

Next, we study problem (ii), the renormalization of the 
relativistic spin-orbit interaction due to the correlation. 
For simplicity, we consider only the z component, 

H,o = (A/2) (£V^) , e^a' = n^^ - ria^ - 7i3t + "3i, (2) 

where V is the 3x3 z component of the angular mo- 
ment matrix of ^ = 1, and a^ is the 2x2 z component 
of Pauli matrix. ^ = 2(3) denotes the pi(p_i) orbital. 
We add Eq.© to Eq.([T|), where the parameters in Eq.([T]) 
are taken as the same values of Figs. [Hb)-(d). For an 
Fe atom, the realistic value of the spin-orbit interaction 
is A = 75 mcV [131. We also show the results with the 
smaller value of A = 40 meV for comparison. Consider- 
ing that the value of A is 2 orders of magnitude less than 
that of impurity energies e j , it is not surprising that the 
temperature dependence of the square of the magnetic 
moment ((Af|)^), the susceptibility Tx^, and the occu- 
pation number (n^) are nearly the same as those of Figs. 
[T]Jb)-(d). In contrast, as displayed in Fig. [^ a nonzero 
spin-orbit correlation function {£^a^) appears when the 
spin-orbit interaction within the t2g orbitals of the Fe im- 
purity is included. At temperature T= 360 K, the lower 
limit of our present calculations, we have {i^a^)= -0.44 
with A = 75 mcV, and {Pa'')^ -0.3 with A = 40 meV. 

To compare the QMC results with the experiment, we 
calculate the spin Hall angle 7s as follows. Since we con- 
sider only two t2g orbitals with the z-component of the or- 
bital angular moment £^ = ±1, the spin-orbit interaction 
within the t2g orbitals gives rise to the difference in the 
occupation numbers between the parallel (np) and an- 
tiparallel (jlap) states of the spin and angular momenta. 
These occupation numbers are related to the phase shifts 
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FIG. 2. (color online) Temperature dependence of spin-orbit 
interactions for an Fe impurity in an Au host. The parameter 
A in Eq. ([2]) is taken as a realistic 75 meV [13, and as 40 
nieV for comparison. 



Sp and Sap, respectively as np(^AP) = Sp(ap)/t^- These 
quantities can be estimated as 

7r(£V^) = Sp - Sap, 7r(n2) + nins) = Sp + Sap, (3) 

which are given in Fig. 2 and 1(d), respectively. Putting 
{i^a')= -0.44 for A= 75 meV, and (712)= (713)= 0.65, we 
obtain Sp= 1.35 and Sap= 2.73. Taking into account 
the estimate 5i ^ 0.1 of the phase shift for p-wave scat- 
tering, and applying the equation [J| 7s = 6/m[(e~^"'i — 
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')] /[25-15cos2(5p-10cos2(5Ap],thespin 
Hall angle is thus obtained as 7s = 0.055, comparable to 
that observed in recent experiment. 

To further understand how correlations renormalize 
the effective relativistic spin-orbit interaction, we will 
now vary the correlation energy and the spin polariza- 
tion, while keeping the temperature T = 360 K and spin- 
orbit interaction A = 75 meV fixed. 

Effect of correlation U: For each U, we keep J/U = 
0.9/5 and U = U' +2J. In Fig. [2]with ((Af|-t-M|)/2) == 
0, we see that the spin-orbit correlation function \{i^a^)\ 
dramatically decreases as U decreases. This demon- 
strates clearly how the correlation U renormalizes the 
relativistic spin-orbit interaction. The results also ex- 
plain past discrepancies in the calculated orbital angular 
momenta for an Fe impurity in an Au host: A large value 
is calculated with U = 5 eV [J| , while much smaller ones 
are obtained with [/ = 3 or 2 eV [1, [11 . The virtual 
bound states of t2g have majority and minority parts, 
whose contributions to {i^a^) arc opposite in sign. In- 
creasing \J will push up the minority part away from the 
Fermi energy, decrease its contribution to {t^a^), and 
hence increase the total {V^u^). 

Effect of spin polarization: To see clearly the effects 
of quantum fluctuations, it is instructive to study the ef- 
fects of polarizing the Fe impurity spin, which can quench 
quantum fluctuations of the spins and orbitals. We apply 
a magnetic field h by hand on the orbitals l(z;^), 2(pi), 
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FIG. 3. (color online) For an Fe impurity in an Au host, 
we show the spin-orbit interaction (l^a^) as a function of 
correlation energy U and spin polarization ({M^ + Ml )/2). 
The QMC results of (^^o"^) are shown as red squares. In the 
phase diagram plane of U to ((Af| + M|)/2), the different 
colors denote the different values of {£^a^). Here, T = 360 K 
and A = 75 meV are fixed. 



and 3(p_i) of the Fc impurity, 

Hf = -h (?iit - nil + "-2t - "24, + "-3t - "-3i) , (4) 

where h may be of a fraction of eV in order to polarize 
the spins at the temperature of 360 K. We add Eq. (|3]) to 
the previous Hamiltonian with the parameters in Eq.([T]) 
taken to be the same as in Fig. [2] It is found that, 
in Fig. [3] with U = 5 eV, when h is set up to 0.8 eV 
to make the spins (M|) of ^ = l(z^), 2(pi) and 3(p_i) 
gradually polarize near the saturated value 1, the spin- 
orbit correlation function K^^cr^)! gradually decreases to 
0. Thus, the quantum fluctuation of the spin is crucially 
important in producing an even larger spin Hall angle 
than the Hall angle of the AHE. What happens is that 
the quantum fluctuation of the spin leads to the motional 
narrowing of the spin-orbit coupled virtual bound states, 
reducing the hybridization energy with the conduction 
electrons and hence the energy width A. So, the Hall 
angle will increase with decreasing A . 

Phase diagram: Combining the above two parameters. 



our QMC results for the spin-orbit correlation function 
{£^a^) are noted as red squares in Fig. [31 In the phase 
diagram plane of the correlation U to the spin polariza- 
tion ((M| -hM|)/2), as also noted in Fig. H the different 
colors note the different values of (I'^a^). It is seen that 
the areas with large spin-orbit correlation should possess 
strong correlation energy U and weak spin polarization. 
The measurement of the orbital angular momentum by 
x-ray MCD involves the application of a sum rule for 
atomic wavefunctions: In order to eliminate interference 
between different impurities, it is necessary to polarize 
with an external magnetic field at low temperatures [7|. 
The spin Hall angle measurement [3| was done in a spin- 
unpolarized state and at room temperature. From our 
phase diagram in Fig. |31 it is suggested that the orbital 
angular momentum in a spin polarized state could be 
much smaller than that in a spin-unpolarized state. We 
emphasize that this phase diagram is obtained for tem- 
peratures comparable to room temperature, as is relevant 
for spin Hall applications; for quantitative comparison of 
very low temperature measurements, such as x-ray MCD 
and phase decoherence, the scales can be renormalized. 

To conclude, we have studied the SHE for an Fe impu- 
rity in Au host by QMC simulation of a multiorbital An- 
derson impurity model. Our phase diagram with respect 
to correlation and spin polarization shows that large spin- 
orbit correlations occur for strong correlation U and weak 
spin polarization, which, in principle, may resolve the 
current, somewhat confusing, situation of experiments 
and calculations. More generally, we can, for the first 
time, quantify an essential difference between the AHE 
and the SHE. While we have focused on the specific case 
of Fe in Au, the results suggest a more general mech- 
anism, and constraints on parameters, for obtaining a 
large SHE by impurities. 
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